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Buoy-mooring  platforms  are  advantageous  for  time-series  validation  and  vicarious 
calibration  of  ocean  color  satellites  because  of  their  high  temporal  resolution  and 
ability  to  perform  under  adverse  weather  conditions.  Bio-optical  data  collected  on 
the  Bermuda  Testbed  Mooring  (BTM)  were  used  for  comparison  with  satellite  ocean 
color  data  in  an  effort  to  further  standardize  sampling  and  data  processing  methods 
for  high  quality  satellite-mooring  comparisons.  Average  percentage  differences  be¬ 
tween  satellite-measured  and  mooring-derived  water  leaving  radiances  were  about 
20%  at  the  blue  wavelengths,  decreasing  to  as  low  as  11%  in  the  blue-green  to  green 
wavebands.  Based  on  a  series  of  data  processing  methods  and  analyses,  recommenda¬ 
tions  concerning  rigor  of  quality  control  for  collected  data,  optimal  averaging  of  high- 
frequency  data,  sensor  self-shading  wind  corrections,  and  instrumentation  placement 
requirements  are  given  for  the  design  and  application  of  optical  moorings  for  ocean 
color  satellite  validation.  Although  buoy-mooring  platforms  are  considered  to  be 
among  the  very  best  methods  to  validate  ocean  color  satellite  measurements,  match¬ 
up  discrepancies  due  to  water  column  variability  and  atmospheric  corrections  re¬ 
main  important  issues. 
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1.  Introduction 

No  other  type  of  observational  platform  has  had  a 
greater  impact  on  large-scale,  synoptic  measurements  of 
the  marine  biological  environment  than  ocean  color  sat¬ 
ellites.  Results  obtained  from  ocean  color  imaging  satel¬ 
lites  such  as  the  Nimbus-7  Coastal  Zone  Color  Scanner 
(CZCS),  the  Sea-viewing  Wide  Field-of-view  Sensor 
(SeaWiFS),  and  the  Moderate-Resolution  Imaging 
Spectroradiometer  (MODIS)  have  greatly  expanded  our 
knowledge  of  global  distributions  of  phytoplankton,  ocean 
primary  productivity,  and  biogeochemical  fluxes  (e.g., 
Aiken  et  al.,  1992;  Hooker  et  al.,  1992;  Behrenfeld  and 
Falkowski,  1997;  Campbell  et  al.,  2002;  Behrenfeld  et 
al.,  2005).  However,  determinations  of  water  column 
properties  from  satellite  platforms  are  complicated  and 
major  improvements  and  increased  information  content 
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are  still  needed  for  many  operational  and  application  pur¬ 
poses.  The  optimal  utilization  of  ocean  color  imagers 
mounted  on  satellite  platforms  is  constrained  by  several 
factors,  including  the  quality  and  interpretation  of  the 
radiometric  data  sets.  Of  utmost  importance  to  data  qual¬ 
ity  are  the  accuracy  of  the  corrections  (including  quanti¬ 
fication  of  the  degradation  of  optical  detectors  and  at¬ 
mospheric  effects)  and  the  efficacy  of  ocean  color  algo¬ 
rithms  that  are  required  to  relate  satellite  ocean  color 
measurements  to  in-water  quantities  (O’Reilly  et  al., 
1998;  Fargion  and  Mueller,  2000;  Hooker  and  McClain, 
2000;  Pinkerton  et  al.,  2003;  lOCCG,  2006).  Thus,  it  is 
critical  that  remote  sensing  platforms  be  complemented 
with  continuous  in  situ  measurements  to  generate  and 
refine  the  necessary  algorithm  input  parameters  for  bio- 
optical  products  (i.e.  Dickey  et  al.,  2006). 

One  of  the  most  successful  approaches  to  validating 
and  vicariously  calibrating  remotely  sensed  ocean  color 
data  has  been  the  use  of  in  situ  optical  measurements,  i.e. 
ship-based,  moored  buoys,  drifters  and  profiling  floats 
(Chavez  et  al.,  1997;  Fargion  and  Mueller,  2000;  Hooker 
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and  McClain,  2000;  Dickey  et  al.,  2001,  2006;  Strutton 
et  al.,  2001;  Clark  et  al.,  2003;  Dickey,  2003;  Kuwahara 
et  al.,  2003;  Pinkerton  et  al.,  2003).  In  particular,  the 
Eulerian  buoy-mooring  observational  platform  has  proven 
to  be  an  effective  vehicle  for  both  validation  and  vicari¬ 
ous  calibration  of  ocean  color  satellites.  Buoy-mooring 
platforms  are  capable  of  sampling  under  even  the  most 
adverse  of  weather  and  sea-state  conditions  with  high  tem¬ 
poral  resolution  for  months  at  a  time  before  the  optical 
sensors  and  system  require  servicing  (Clark  et  al.,  1997; 
Dickey  et  al.,  2001, 2006).  On  the  other  hand,  ocean  color 
satellites  are  effectively  restricted  during  even  the  slight¬ 
est  cloud  coverage  and  wind  generated  white-cap  condi¬ 
tions.  Multi-year  deployments  with  periodic  recoveries 
and  re-deployments  of  buoy-mooring  systems  fitted  with 
optical  instrumentation  (hereafter  referred  to  as  optical 
moorings)  have  become  possible  and  successful  as  a  re¬ 
sult  of  technological  advances  in  hardware,  power 
sources,  data  storage,  telemetry  capabilities,  and  copper 
anti-fouling  devices  (Chavez  et  al.,  2000;  Manov  et  al., 
2004). 

Optical  moorings  have  formed  the  foundation  of  sev¬ 
eral  long-term  ocean  monitoring  projects  and  validation 
and/or  vicarious  calibration  efforts  (Chavez  et  al.,  1997; 
Clark  et  al.,  1997,  2003;  Ishizaka  et  al.,  1997;  Kishino  et 
al.,  1997;  Dickey  et  al.,  1998,  2001,  2006;  Letelier  et  al., 
2000;  Barnes  et  al.,  2001 ;  Strutton  et  al.,  2001 ;  Pinkerton 
et  al.,  2003;  Spada  et  al.,  2007).  Moorings-of-opportu- 
nity  can  provide  comprehensive,  interdisciplinary  in  situ 
observations  complementing  satellite  remotely  sensed 
data  for  deriving  synergistic  descriptions  of  oceano¬ 
graphic  features  and  their  evolution  in  space  and  time,  as 
well  as  providing  optical  validation  and,  potentially,  vi¬ 
carious  calibration  data.  Optical  moorings  have  helped 
facilitate  ocean  color  satellite  validation  efforts  by  pro¬ 
viding  high  frequency,  ground-truth  measurements  in  a 
variety  of  open  ocean  locations;  they  have  already  col¬ 
lectively  influenced  correction  techniques  and  algorithm 
development. 

Although  the  potential  of  optical  moorings  has  been 
demonstrated,  intercomparisons  between  different  meas¬ 
urement  and  methodological  techniques  have  yet  to  be 
standardized.  For  example,  protocols  to  calculate  relevant 
in  situ  water-leaving  radiance,  L^(A),  and  normalized 
water-leaving  radiance,  L„„(A),  have  been  recommended, 
but  detailed  descriptions  of  instrumentation  placement, 
sampling  frequency,  site-specific  limitations,  etc.  have  yet 
to  be  standardized  and  are  thus  inconsistent  between  stud¬ 
ies  (Kuwahara  et  al.,  2003).  In  order  to  assure  that  radio- 
metric  and  bio-optical  data  acquired  from  optical  moor¬ 
ings  meet  standards  of  quality  and  accuracy,  clear  and 
well-documented  sampling  and  data  processing  methods 
need  to  be  standardized. 

The  purpose  of  this  study  is  to  make  recommenda¬ 


tions  for  the  design  and  application  of  optical  moorings- 
of-opportunity  for  ocean  color  satellite  ground-truthing 
purposes  on  a  global  scale.  Specifically,  we:  (1)  demon¬ 
strate  that  moorings-of-opportunity  are  viable  platforms 
for  ocean  color  satellite  validation  by  comparing  in  situ 
optical  mooring  and  satellite  ocean  color  match-ups  of 
L^(A),  (2)  discuss  the  quality  control  and  data  analysis 
protocols  necessary  for  quality  measurements,  and 
(3)  examine  different  ways  to  optimally  deploy  instru¬ 
mentation  to  minimize  errors. 

2.  Materials  and  Procedures 

Data  relevant  to  this  study  were  collected  on  the 
BTM,  which  was  designed  and  configured  to  provide  high 
frequency  in  situ  measurements  at  a  fixed  location  in  deep 
waters  off  Bermuda  for  oceanographic  science  (Dickey 
et  al.,  2001).  The  BTM  is  located  at  roughly  31°N,  64°W, 
about  80  km  southeast  of  Bermuda  in  waters  of  -4567  m 
depth.  BTM  instrumentation  of  primary  interest  to  the 
present  study  includes  a  chlorophyll  fluorometer  at  35  m 
and  radiometers  for  upwelling  radiance  and  downwelling 
irradiance.  Each  of  the  latter  systems,  co-located  at  14 
and  21m  depths,  utilize  Satlantic,  Inc.  OCT200  and  OCR- 
200  radiometers  that  measure  downwelling  irradiance, 
Ejj(A,  z),  and  nadir  upwelling  radiance,  L^(A,  z),  at  seven 
wavelengths  centered  at  412,  443,  490,  510,  555,  665, 
and  683  nm  (spectral  bandwidth  is  10  nm).  These  wave¬ 
lengths  are  compatible  with  those  of  the  SeaWiFS  ocean 
color  imager  plus  683  nm.  In  addition  to  the  subsurface 
radiometer  packages,  an  OCI-200  was  mounted  on  the 
BTM  surface  buoy  just  above  the  sea  surface  for 
Ejj(A,  0"^)  and  a  Satlantic,  Inc.  Stordat-7  radiometer  was 
deployed  at  5  m  depth  to  measure  5m)  at  the  same 
seven  wavelengths  as  those  listed  above.  Sampling  by  the 
OCI/OCR-200  radiometers  was  done  for  45  seconds  at  6 
Hz  every  hour.  The  sampling  rate  for  the  Stordat-7  was  2 
Hz  every  hour  for  30  seconds.  All  sensors  were  factory 
calibrated  following  the  standards  set  down  by  the  ocean 
color  community  (Mueller  and  Austin,  2003).  Although 
pre-  and  post-calibrations  of  all  radiometers  were  con¬ 
ducted,  we  utilized  the  post-calibration  for  the  present 
study  since  calibration  drift  was  minimal  for  all  instru¬ 
ments  at  all  wavelengths  (<1.5%).  The  data  presented  here 
were  collected  during  the  12th  deployment  of  the  BTM, 
between  30  July  and  5  November  1999. 

2.1  Data  quality  assurance/quality  control  (QA/QC) 

A  data  pre-processing  routine  was  developed  to  se¬ 
lect  BTM  data  deemed  appropriate  for  analysis.  This  rou¬ 
tine  involved  calculation  of  the  calibrated  measurements 
for  the  entire  deployment  at  each  wavelength  for  each 
instrument  at  each  depth  and  a  series  of  quality  assur¬ 
ance  and  quality  control  (QA/QC)  analyses,  described 
here  (Fig.  lA). 
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C.  Error  Analyses 


Fig.  1.  Flow-chart  showing  the  procedures  used  during  this  study:  data  quality  assurance/quality  control  (QA/QC),  post-process¬ 
ing,  and  error  analyses. 


(1)  The  time-series  data  were  inspected  and  scruti¬ 
nized  for  obvious  data  dropouts  from  instrument/power 
failure  and  symptoms  of  salt  or  other  depositions  on  the 
optical  windows  of  the  radiometers.  Subsurface  radiomet¬ 
ric  data  were  subjected  to  the  test  for  fouling  following 
work  by  Abbott  and  Letelier  (1998),  who  suggested  that 
data  are  affected  by  bio-fouling  when  the  upwelling  radi¬ 
ance  ratio  of  683/555  nm  exceeds  0.1,  which  is  a  value 
appropriate  for  clear  oligotrophic  water  masses  such  as 
those  found  near  Bermuda.  In  the  present  study,  ratio  val¬ 
ues  for  each  depth  (5,  15,  and  21  m)  were  well  below  the 
0.1  threshold  level,  suggesting  negligible  influence  from 
bio-fouling  and  effectiveness  of  precautionary  anti-foul- 
ing  methods  (i.e.  copper  shutters;  Manov  et  al.,  2004). 
No  data  points  showed  signs  of  obvious  instrument  deg¬ 
radation  or  failure. 

(2)  The  time-series  of  normalized  spectra  and  wave- 
length-to-wavelength  ratios  from  each  surface  and  sub¬ 
surface  data  set  were  analyzed  to  determine  the  consist¬ 
ency  and  shape  of  each  spectrum.  The  purpose  of  this 
exercise  is  to  determine  whether  the  wavelength-to-wave- 
length  ratio  maintains  consistency  or  “rank-order” 
throughout  the  deployment  as  defined  by  the  relative 
magnitude  of  normalized  irradiance  at  all  wavelengths, 
and  to  determine  whether  the  wavelength-to-wavelength 
ratio  values  fall  within  limits  defined  by  clear-sky  mod¬ 
els  of  incident  daylight  (Kuwahara  et  al.,  2003).  It  is  ex¬ 
pected  that  the  collected  data  should  maintain  a  consist¬ 
ent  normalized  rank-order  unless  the  instrument  sensor 


is  damaged  or  improperly  calibrated.  The  results  obtained 
from  the  time  series  of  interest  indicate  that  consistent 
wavelength-to-wavelength  order  was  maintained;  this 
implies  that  sensor  diode  integrity  and  reliable  pre-  and 
post-deployment  calibration  prevailed. 

(3)  For  incident  spectral  irradiance  above  the  sea 

surface,  0"^),  a  comparison  of  the  magnitudes  of 

measured  E(](A,  0"^)  with  clear-sky  model  clear-sky) 
estimates  by  Frouin  et  al.  (1989)  were  calculated  for  the 
solar  zenith  angle  at  each  measurement  time  for  the  en¬ 
tire  deployment.  Data  were  labeled  as  suspect  if  the  inci¬ 
dent  measurements  exceeded  a  threshold  factor  of  1.25  x 
E(](A,  clear-sky).  The  factor  1.25  allows  measured  spec¬ 
tral  irradiances  to  moderately  exceed  calculated  clear-sky 
irradiances  due  to  scattering  from  clouds.  None  of  the 
daily  average  EJ^X,  0"^)  values  failed  this  clear-sky  qual¬ 
ity  control  test. 

(4)  Self-shading  by  the  radiometers  was  corrected 
following  the  methods  presented  by  Gordon  and  Ding 
(1992).  The  error  associated  with  radiometer  self-shad¬ 
ing,  s{X),  can  be  represented  as: 

e(X)  =  [Ej(X)  -  L  “(A)]/L„T(A),  (1) 

and 

e{X)  =  1  -  exp[-k  aj(A)r],  (2) 

where  Ly^(A)  is  radiance  corrected  for  self-shading  and 
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Lu’^(A)  is  uncorrected  radiance,  a^{X)  is  the  total  absorp¬ 
tion  coefficient,  r  is  the  radius  of  the  instrument  housing 
(r  =  0.44  m),  and  k  =  2/tan Oqw  (^Ow  the  refracted  solar 
zenith  angle).  This  method  assumes  that  absorption  domi¬ 
nates  over  scattering  processes,  which  is  true  for  Sargasso 
Sea  waters.  The  absorption  coefficient  was  estimated  us¬ 
ing  the  model  presented  by  Morel  (1991): 

a, (A)  =  [a„(A)  +  0.06a,u*(A)ChlO®5][l  +  0.2y(A)],  (3a) 
where 

y(A)  =  exp[-0.014(A  -  440)],  (3b) 

a^(A)  is  the  absorption  coefficient  of  pure  water  (Pope 
and  Fry,  1997),  aj,jji*(A)  is  the  chlorophyll-specific  ab¬ 
sorption  coefficient  (Prieur  and  Sathyendranath,  1981), 
and  Chi  is  the  BTM-measured  chlorophyll-a  concentra¬ 
tion.  Because  values  of  a^{X)  are  extremely  low  in  the 
Sargasso  Sea,  radiometer  self-shading  effects  were  small, 
with  £(X)  ~  1%.  We  assumed  that  buoy-shading  effects 
were  negligible;  radiometers  were  mounted  well  below 
the  base  of  the  buoy  (z  >  5  m)  and  on  an  arm  extending 
away  from  the  centerline  of  the  mooring.  In  addition,  sat¬ 
ellite  and  mooring  match-up  times  used  in  this  study  were 
all  within  two  hours  of  Bermuda  local  noon,  reducing  the 
effects  of  the  buoy  shading. 

2.2  Data  post-processing 

After  performing  the  QA/QC  protocols  for  moored 
radiometric  data,  SeaWiFS  and  mooring  data  match-up 
times  were  identified.  Successful  match-ups  are  most 
dependent  on  clear  sky  conditions  and  equivalent  spatial 
coverage  (i.e.  mooring  occurring  in  a  pixel).  Of  the  en¬ 
tire  101-day  BTM  radiometric  time  series  between  30  July 
and  5  November  1999,  23  suitable  SeaWiFS  satellite 
match-ups  were  available,  or  roughly  23%  of  collected 
data.  For  comparison,  PlyMBODy  optical  data,  collected 
in  a  relatively  cloudy  oceanic  region  (English  Channel), 
matched  up  less  than  6%  of  the  time  during  a  ~10  month 
deployment  and  about  10%  of  MOBY  data  are  realized 
for  SeaWiFS  match-ups  (Pinkerton  et  al.,  2003).  These 
differences  underscore  both  the  limitations  and  benefits 
of  using  optical  moorings-of-opportunity  for  satellite  vali¬ 
dation  and  vicarious  calibration,  as  moorings  provide 
autonomously  collected,  long,  continuous  time  series  from 
which  satellite  match-up  data  may  be  selected. 

For  the  purpose  of  ocean  color  satellite  validation  it 
is  necessary  to  average  individual  high-frequency  meas¬ 
urements.  The  averaging  of  high  frequency  data  accounts 
for  fluctuations  in  the  incident  light  field,  such  as  those 
caused  naturally  by  wave-focusing  (e.g.,  Zaneveld  et  al., 
2001;  Zheng  et  al.,  2002),  and  allows  comparisons  to  be 
made  with  satellites  using  a  single,  averaged  measure¬ 


ment  for  a  particular  satellite  image.  For  the  BTM  data 
sets,  three  different  averaging  techniques  were  tested.  For 
each  wavelength  at  each  depth:  (1)  daily  averages  were 
calculated  using  high-frequency  data  from  0700  to  1900 
Atlantic  Standard  Time  (AST),  (2)  averages  were  com¬ 
puted  for  daily  noon  between  1000  and  1400  AST,  and 
(3)  1-hr  averages  were  computed,  centered  around  the 
time  nearest  to  the  SeaWiFS  satellite  overpass  times.  Sat¬ 
ellite  overpass  times  were  all  within  1-hr  of  mooring  data 
collection. 

The  vertical  diffuse  attenuation  coefficients  for 
downwelling  irradiance  and  upwelling  radiance,  Kjj(Z,  z) 
and  Kl(A,  z),  respectively,  were  determined  directly  from 
QA/QC’ed  radiometric  measurements.  Given 
downwelling  irradiance  measured  by  radiometers  at  two 
different  depths,  the  diffuse  attenuation  coefficient  for 
downwelling  irradiance  at  the  midpoint  of  the  two  depths 
is  given  by: 


(4a) 

1  j^Ed(Z,Z2) 

Az  Ejj(/l,Zi)’ 

(4b) 

where  depths  Z2  >  Zj.  K(](Z,  18m)  was  computed  using 
measured  E(](Z,  z)  at  15  and  21  m.  Eurther  QA/QC  analy¬ 
ses  are  performed  on  calculated  z);  K^{X,  z)  should 
be  greater  than  the  absorption  coefficient  of  pure  water 
at  each  wavelength  (Pope  and  Pry,  1997).  Throughout  the 
time  series,  K^{X,  18m)  values  for  X  =  665  and  683  nm 
failed  the  pure  water  absorption  coefficient  test,  likely 
due  to  the  low  signal-to-noise  ratio  for  Ejj(Z,  z)  at  the  red 
wavelengths.  Similar  to  KJ^X,  z),  given  two  underwater 
upwelling  radiance  sensors,  the  vertical  attenuation  co¬ 
efficient  for  upwelling  radiance  at  the  center  of  that  depth 
interval  is  calculated  by: 


-£[1"Lu(A,z)], 

(5a) 

Az  L^(Z,Zi)' 

(5b) 

Kl(A,  10m)  was  computed  from  measurements  of 
Lu(Z,  5m)  and  L^^{X,  15m),  whereas  Kl(A,  18m)  was  ob¬ 
tained  using  LJ^X,  z)  measurements  at  15  and  21m,  and 
finally  Kl(A,  13m)  was  calculated  from  L^(Z,  5m)  and 

L^(X,  21m). 

Water  leaving  radiance,  L„(Z),  was  determined  by 
extrapolating  upwelling  radiance  measured  at  depth  z  to 
the  surface  following: 
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Table  1.  Comparisons  of  0“)  calculated  using  Eq.  (6a)  and  vertical  diffuse  attenuation  coefficients  for  upwelling  radiance 
and  downwelling  irradiance  computed  at  and  extrapolated  from  various  depths.  Each  table  cell  shows  results  at  three  wave¬ 
lengths:  412,  490,  and  555  nm,  respectively  computed  using  averages  between  0700  and  1900  AST.  Percentage  differences 
were  calculated  using  Eq.  (7)  and  averaged  over  the  time  series.  Differences  of  less  than  10%  are  indicated  in  bold. 


LJAO  )  extrapolated  using: 

L,(A,0  )  extrapolated  using: 

K,(18m)L„(5m) 

K,(18m)  LT15m) 

K,(18m)  LT21m) 

KT10m)L„(5m) 

KTlSm)  LT15m) 

KL(13m)  L,(21m) 

K,(18m)  L,(5m) 

6.7% 

15% 

-3.3% 

-13% 

-4.6% 

— 

-6.3% 

0.16% 

3.2% 

-23% 

-0.05% 

-37% 

-32% 

19% 

-51% 

10% 

K,(18m)  LjlSm) 

-6.7% 

7.9% 

-10% 

-20% 

-11% 

6.3% 

— 

6.5% 

9.5% 

-16% 

6.3% 

37% 

5.7% 

55% 

-14% 

47% 

K,(18m)  L,(21m) 

-15% 

-7.9% 

-18% 

-28% 

-19% 

-0.16% 

-6.5% 

— 

3.0% 

-23% 

-0.21% 

32% 

-5.7% 

50% 

-20% 

41% 

K^dOm)  L,(5m) 

3.3% 

10% 

18% 

-9.8% 

-1.2% 

-3.2% 

-9.5% 

-3.0% 

— 

-26% 

-3.2% 

-19% 

-55% 

-50% 

-57% 

-8.8% 

KL(18m)  LJlSrn) 

13% 

20% 

28% 

9.8% 

8.6% 

23% 

16% 

23% 

26% 

— 

22% 

51% 

14% 

20% 

57% 

60% 

KTlBm)  L,(21m) 

4.6% 

11% 

19% 

1.2% 

-8.6% 

0.05% 

-6.3% 

0.21% 

3.2% 

-22% 

— 

-10% 

-47% 

-41% 

8.8% 

-60% 

L,(A,0-)=L,(A,z)e'^(^'^)^  (6a) 

and  through  the  surface  according  to: 

L„(A)  =  0.543L,(A,  01,  (6b) 

where  only  nadir-viewing  geometry  is  considered  (the 
surface  refiectance  becomes  independent  of  wind  speed) 
and  the  constant  vaiue,  0.543,  corresponds  to  theupward 
transmittance  term  (Austin,  1974).  LJA,  01  and  L„(A) 
werecomputed  12  different  ways  using  Eqs.  (6a)  and  (6b): 
(l)-(3)  Kl(A,  10m)  and  L^A,  z),  (4)-(6)  Kl(A,  13m)  and 
L^A,  z),  and  (7)-(9)  Kl(A,  18m)  and  LJA,  z),  where z  = 
5,  15,  and  21  m;  K(j(A,  18m)  was  substituted  for  Kl(A,  z) 
and  (10)  L^A,  5m),  (11)  L^A,  15m),  and  (12)  L^A,  21m). 
For  this  study  we  present  L„(A)  rather  than  normaiized 
water-ieaving  radiance,  in  order  to  evaiuate  ef¬ 

fects  of  soiar  zenith  angie  variabiiity. 

2.3  SeaWiFS  data  processing 

SeaWiFS  data  were  provided  by  theNASA  SIM  BIOS 
(National  Aeronautics  and  SpaceAdministration  Sensor 
Intercomparison  for  M  arine  B iological  and  Interdiscipli¬ 


nary  Ocean  Studies)  projectfor  a  5  x  5  pixel  box  (roughly 
25  sq.  km)  centered  on  the  nominal  BTM  location.  The 
5x5  pixel  box  was  selected  as  opposed  to  a  1  x  1  pixel 
box  to  maximize  satellite  mach-ups  due  to  environmen¬ 
tal  variability  at  the  BTM  site  and  to  provide  more  robust 
data.  A  total  of  1,251  SeaWiFS  files  (global  and  site-spe¬ 
cific  local  area  coverage)  were  processed:  246  of  these 
files  passed  to  Level  2  SeaWiFS  products  and  were  con¬ 
sidered  useable.  After  reducing  data  points  to  those  that 
had  at  least  one  valid  pixel  in  the  5  x  5  box,  87  match-ups 
remained.  Of  these  87  total  match-ups,  46  local  area  cov¬ 
erage  data  poi  nts  remai  ned  once  the  redundant  data  from 
gl obal  area  coverage  were  removed.  Twenty-three  of  these 
data  points  matched-up  with  the  BTM  data  collection 
period  and  were  used  for  the  comparative  analysis  with 
the  BTM  radiometer  data.  Due  to  frequent  failure  of  the 
pure  water  absorption  coefficient  test,  the  red  wavelengths 
of  665  and  683  nm  were  not  used  in  the  analyses  pre¬ 
sented  here. 

3.  Assessment  and  Results 

In  the  present  study,  downwelling  irradiance  and 
upwelling  radiance  were  measured  on  a  mooring  at  mul¬ 
tiple  depths,  allowing  us  to  establish  the  requirements  for 
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Fig.  2.  Mean  spectra  averaged  over  the  time  series  (n  =  23;  aii  spectra  shown  in  grey)  of:  (a)  BTM -measured  15m); 
(b)  Kj(A,  18m)  computed  using  Eqs.  (4a)  and  (4b);  (c)  LJA,  0“)  extrapoiated  using  Kj(A,  18m)  and  LJA,  z)  (Eq.  (6a)); 
(d)  BTM  -measured  LJA,  15m);  (e)  Kl(A,  z)  computed  using  Eqs.  (5a)  and  (5b);  (f)  L,j(A,  0“)  extrapoiated  using  Kl(A,  Zj)  and 
LJA,  Zj)  (Eq.  (6a)).  Comparisons  between  0“)  extrapoiated  using  (g)  Kc|(A,  18m)  and  L,j(A,  z^)  vs.  LJA,  Zj);  (h)  Kl(A,  Zj) 
versus  Kl(A,  Zj)  and  L,j(A,  z)  as  in  Tabie  1;  (i)  Kc|(A,  18m)  and  Lu(A,  z)  versus  Kl(A,  13m)  and  L,j(A,  21m).  Waveiength 
notations  in  figure  are  suppressed. 


reliable  measurements  necessary  for  validation  of  ocean 
color  satellite  data. 

3.1  Vertical  diffuse  attenuation  coefficients,  KJ f,  z)  and 

Kl(1  z) 

The  vertical  diffuse  attenuation  coefficients  for 
upwelling  radianceand  downwelling  irradiance,  Kl(A,  z) 
and  K(j(A,  z),  evaluated  for  multiple  depths  were  used  to 
extrapolate  L^j(A,  Oi  (Fig.  3).  Comparisons  were  made 
between  derived  values  of  L^j(A,  O^j  (Table  1)  to  deter¬ 
mine  the  optimal  placement  of  radiometers  for  satellite 
validation  purposes.  Percentage  differences  were  com¬ 
puted  using  the  following  equation: 


%diff(A)  =  100x 


0.: 

5[l,(a,o-)^-l4a,o-)^ 

where  the  subscripts  1  and  2  specify  two  different  meth¬ 
ods  of  computing  LJA,  Oi,  e.g.,  1:  using  Kl(A,  10m)  and 
L  JA,  5m)  and  2:  using  K(j(A,  18m)  and  L  JA,  5m).  Note 
that  these  comparisons  involved  mooring  measurements 
only.  In  general,  derivations  of  LJA,  0^)  using  K(j(2„  18m) 


with  L  JA,  5m),  L  JA,  15m),  and  LJA,  21m)  (Eq.  (6a)) 
were  comparable  with  average  percentage  differences  of 
near  10%  for  all  wavelengths  except  555  nm  (Table  1; 
Fig.  2).  As  expected,  larger  deviations  were  found  when 
comparing  vertical  diffuse  attenuation  coefficients  using 
5  m  extrapolated  to  21  m.  Also  as  expected,  for  deriva¬ 
tions  of  L^j(A,  Oi  using  Kl(A,  z)  with  L^j(A,  5m;  15m;  21m) 
at  different  depths,  differences  were  less  than  1%  when 
comparing  extrapolations  from  the  same  depths  used  in 
the  computation  of  K  l(A,  z)  and  on  the  order  of  10%  for 
extrapolations  from  the  other  two  depths  (data  not  shown). 
For  example,  L^j(A,  Oi  computed  using  Kl(A,  10m)  and 
LJA,  5m)  was  comparable  to  L^j(A,  oi  computed  using 
Kl(2,,  10m)  and  L  JA,  15m)  because  Kl(A,  10m)  was  de¬ 
rived  from  measurements  of  L  JA,  5m)  and  L  JA,  15m). 

Comparisons  between  L^j(A,  0^)  derived  using 
Kl(A,  z),  wherez  =  10,  13,  or  18  m,  varied  depending  on 
the  extrapolation  depth.  Lowest  percentage  differences 
were  found  for  L  JA,  Oi  derived  using  Kl(2,,  10m)  and 
Kl(A,  13m),  i.e.  using  the  two  uppermost  radiometers  and 
the  upper  and  lowest  mounted  radiometer  (Table  1).  The 
largest  deviations  were  found  for  computations  using 
KlU,  18m),  illustrating  that  optical  properties  between 
15  and  21  m  were  differentfrom  those  properties  between 
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Table  2.  Comparisons  between  SeaWiFS-measured  and  mooring-derived  values  of  L„(A).  Each  table  cell  shows  three  different 
averaging  schemes:  daily  averaged  (0700-1900  AST),  4-hr  averaged  (1000-1400  AST),  and  1-hr  averaged  (nearest  satellite 
overpass  time)  differences,  respectively.  Percentage  differences  were  computed  using  Eg.  (8)  and  averaged  over  the  time 
series. 


SeaWiFS-measured:  L„(A)  derived  using  mooring: 


Kj(18m)  L„(5m) 

Kd(18m)  L„(15m) 

Kj(18m)  L,(21m) 

K^flOm)  L  j5m) 

KL(18m)  L„(15m) 

KL(13m)  L  j21m) 

L»(412  nm) 

-56% 

-59% 

-62% 

-54% 

-49% 

-54% 

-30% 

-37% 

-42% 

-25% 

-22% 

-25% 

-23% 

-34% 

-40% 

-17% 

-17% 

-17% 

L„(443  nm) 

-60% 

-61% 

-65% 

-60% 

-50% 

-59% 

-37% 

-41% 

-47% 

-35% 

-24% 

-34% 

-32% 

-39% 

-45% 

-27% 

-20% 

-27% 

L„(490  nm) 

-57% 

-54% 

-57% 

-58% 

-46% 

-57% 

-32% 

-31% 

-35% 

-32% 

-17% 

-30% 

-26% 

-28% 

-34% 

-24% 

-12% 

-23% 

L„(510  nm) 

-64% 

-58% 

-52% 

-66% 

-45% 

-64% 

-42% 

-36% 

-43% 

-45% 

-16% 

-42% 

-38% 

-34% 

-41% 

-38% 

-11% 

-36% 

L„(555  nm) 

-61% 

-43% 

-47% 

-68% 

-34% 

-65% 

-38% 

-13% 

-18% 

-47% 

3.8% 

-42% 

-33% 

-9.6% 

-16% 

-41% 

11% 

-37% 

5  and  15  m.  L^j(A,  0^)  derived  using  K(j(A,  z)  versus  those 
using  Kl(A,  z)  were  comparabie,  except  for  Kl(A,  18m) 
(Table  1;  Fig.  2).  In  general,  differences  computed  at  555 
nm  were  greater  than  those  at  the  other  wavelengths.  Ra¬ 
diance  spectra  were  at  a  minimum  at  this  wavelength;  low 
signal-to-noise  ratios  in  clear  Sargasso  Sea  waters  may 
have  been  the  reason  for  poor  correlations  at  555  nm  (Fig. 
2).  Daily  averaged  (over  daylight  hours),  4-hr  averaged 
(centered  at  noon),  and  1-hr  averaged  (centered  at  over¬ 
pass)  results  were  quite  similar  to  each  other. 

BTM  data  showed  that  K(j(2,,  z)  was  reliable  in  propa¬ 
gating  L  JA,  z)  measurements  at  depth  to  just  below  the 
surface.  Because  K(j(A,  z)  undergoes  a  quality  control 
screening  (pure  water  absorption  coefficient  comparison) 
and  it  is  an  important  variable  in  many  bio-optical  mod¬ 
els  (e.g.,  M  Orel  and  M  aritorena,  2001),  it  may  be  useful 
in  some  operations  to  utilize  K(j(A,  z)  to  extrapolate 
L^j(2„  z)  to  the  surface  only,  when  Kl(A,  z)  is  not  readily 
available.  In  other  words,  if  a  limited  number  of  radiom¬ 
eters  can  be  deployed  or  if  some  sensors  fail,  it  may  be 
desirable  to  utilize  one  L  JA,  z)  sensor  and  two  E(j(iL,  z) 
sensors.  Regardless  of  whether  K(j(A,  z)  orKL(A,  z)  is  (or 
both  are)  used,  sensors  deployed  at  multiple  depths  are 
strongly  recommended  to  assure  against  data  loss  in  the 
event  of  sensor  failure. 


3.2  Water-leaving  radiance 

Water-leaving  radiances  were  computed  using  all 
derived  L  JA,  Oi  quantities  and  Eq.  (6b)  (12  different 
computations).  These  mooring-derived  values  were  then 
compared  to  SeaWiFS-measured  L„(A)  (Fig.  3;  Table  2). 
Percentage  differences  were  computed  according  to  the 
following  equation: 

%diff(A)  =  100  X  [(L„(A)2  -  L„(A)i)/L„(A)i],  (8) 

where  L„(A)2  and  L„(2,)i  are  water-leaving  radiances 
measured  by  SeaWiFS  and  derived  from  BTM  radiomet¬ 
ric  measurements,  respectively.  Comparisons  between 
satellite-and  mooring-derived  L„(A)  were  wavelength 
dependent.  For  the  uppermost  water  column  (5-15  m), 
percentage  differences  generally  increased  with  wave¬ 
length.  Interestingly,  the  opposite  was  true  for  L„(A)  de¬ 
rived  from  LJA,  21m).  The  reason  for  this  is  unknown. 
Pinkerton  et  al.  (2003)  reported  higher  absolute  percent¬ 
age  differences  between  SeaWiFS-measured  and 
PlyM  BODy-derived  L^JA)  at  the  blue  (412-443  nm)  and 
red  (670  nm)  wavelengths  (56  and  77%)  compared  to  the 
blue-green  to  green  wavebands  (490-555  nm  =  -20%), 
likely  due  to  increased  turbidity  in  the  region. 

BTM  results  showed  that  mooring  data  averaged  over 
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Fig.  3.  Time-averaged  water-ieaving  radiance,  L^(A),  spectra 
(n  =  23)  measured  by  SeaWiFS  and  derived  (a)  using 
Kc|(A,  18m)  and  z);  (b)  using  Kl(A,  z^)  and  Zj). 
Comparisons  between  L„(A)  (c)  derived  using  thethreedif- 
ferent  depths  in  (a);  and  (d)  derived  using  the  three  differ¬ 
ent  depth  combinations  in  (b). 


daylight  hours  (0700-1900  AST)  did  not  compare  well 
with  ocean  color  satellite  measurements,  which  are  col¬ 
lected  at  one  instant  in  time.  Absolute  differences  were 
of  the  order  of  50-70%  for  all  wavelengths  at  all  depths 
using  K  l(A,  z)  or  K(j(A,  18m).  On  the  other  hand,  results 
for  averaging  performed  over  4-hrs  and  1-hr  were  im¬ 
proved  as  indicated  by  absolute  differences  of  roughly 
25%  (30-40%)  and  <20%  (20-40%),  respectively  for 
L„(A)  derived  using  Kl(412  nm,  z)  [K(j(412  nm,  18  m)] 
(Fig.  3;  Table  2).  In  all  three  time  averaging  cases,  deri¬ 
vations  using  Kl(A,  18m)  with  L  JA,  15  and  21m)  were 
better  than  those  for  Kl(A,  10m),  Kl(A,  13m),  or 
K(j(A,  18m).  Recall  that  results  from  vertical  diffuse  at¬ 
tenuation  coefficient  comparisons  showed  that  optical 
properties  between  15  and  21  m  were  different  from  those 
properties  between  5  and  15  m. 

4.  Discussion 

In  general,  the  differences  between  satellite- mea¬ 
sured  and  BTM  -derived  L„(2,)  were  comparable  to  those 
found  using  other  optical  moorings,  except  when  aver¬ 
ages  between  0700-1900  AST  were  used.  Average  dif¬ 
ferences  over  the  SeaWiFS  wavebands  were  37%  for 
PlyM  BODy  results  (Pinkerton  et  al,  2003).  Water-leav¬ 
ing  radiance  derived  using  Kl(2„  18m)  and  L^j(A,  15m) 
were  greatly  improved,  with  differences  of  20%  in  the 
blue,  decreasing  to  11%  in  the  green  wavebands.  Various 
error  analyses  associated  with  in  situ  measurements  were 
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Fig.  4.  KjfA,  18m)  and  Kl(2„  18m)  error  analysis  results, 
(a)  Kj(A,  18m)  averaged  over  the  time  series  (n  =  23;  black 
line  with  closed  circles)  plotted  with  Kj(A,  18m)  + 
Kd(A,  18m)  X  20%,  40%,  60%,  80%,  and  100%.  (b)  As  in 
(a)  but  for  Kl (A,  18m).  (c)  and  (d)  Derived  L„(A)  using  ver¬ 
tical  diffuse  attenuation  coefficients  in  (a)  and  (b),  respec¬ 
tively,  and  Eqs.  (6a)  and  (6b).  (e)  and  (f)  Percent  differ¬ 
ences  in  L„(A)  (Eg.  (8)). 


conducted  to  evaluate  the  impact  of:  instrumentation  in¬ 
accuracies  (e.g.,  radiometric  drift,  miscalibrations, 
biofouling,  wavelength  shifts),  mooring  motion,  and  en¬ 
vironmental  effects  (solar  angle  and  wind  speed)  on  the 
discrepancies  related  to  methods  of  estimating  L„(A)  (Fig. 
1C).  Radiance  is  an  extremely  difficult  measurement  to 
make.  For  example,  Chang  etal.  (2003)  found  differences 
of  -20%  between  three  different  in  situ  methods  of  esti¬ 
mating  L„(A). 

4.1  Instrumentation  effects 

Propagation  of  errorsfrom  inaccurate  measurements 
of  E(j(2„  z)  and  L  JA,  z)  caused  by,  e.g.,  radiometer  drift, 
miscalibration,  or  biofouling,  were  investigated  by  com¬ 
puting  K(j(A,  18m)  and  Kl(X  18m)  from  time-averaged 
BTM  measurements  of  E(j(A,  z)  and  L^j(A,  z),  where  z  = 
15  and  21,  then  varying  these  values  of  K(j(2„  18m)  and 
Kl(A,  18m)  from  0  to  +100%  by  steps  of  +20%.  Water¬ 
leaving  radiance  was  derived  from  variable  values  of 
K(j(A,  18m)  and  Kl(A,  18m)  and  BTM -measured 
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Fig.  5.  Results  from  mooring  heave  error  analysis.  L„(A)  spectra  extrapolated  using  averaged  (n  =  23)  (a)  18m  +  Az)  and 

Lu(A,  5m  +Az);  (b)  Kl(A,  10m  +Az)  and  L,j(A,  5m  +Azj;  (c)  Kl(A,  13m  +Az)  and  LyfA,  21m  +Az);  and  (d)  Kl(A,  18m  +Az) 
and  L,j(A,  15m),  and  Az  =  0.25  to  3  m  by  increments  of  0.25  m.  Percentage  differences  in  L„(A)  between  non-depth-shifted 
(Az  =  0  m)  and  depth-shifted  parameters  (Eg.  (8))  (e)-(h)  derived  using  (a)-(d). 


L^j(A,  15m)  (Eqs.  (6a)  and  (6b)).  The  results  of  thisanaly- 
sis  are  shown  in  Fig.  4.  Naturally,  higher  values  of  dif¬ 
fuse  attenuation  coefficients  resulted  In  greater  changes, 
therefore  the  computation  of  L„(A)  using  Kl(A,  18m)  was 
much  more  sensitive  to  changes  than  those  using 
K(j(A,  18m),  and  longer  wavelengths  were  most  affected. 
Moreover,  Increases  as  opposed  to  decreases  In 
K(j(A,  18m)  and  Kl(A,  18m)  resulted  In  higher  percent¬ 
age  differences  between  derived  L„(A).  Generally,  a  20% 
errorin  the  computation  of  K(j(A,  z)  and  Kl(A,  z)  resulted 
In  <10%  errors  In  the  derivation  of  L„(A),  however  up  to 
60%  errors  were  found  when  Kl(A,  z)  was  100%  of  Its 
true  value  (blue  wavelengths).  These  errors  Increase 
exponentially  with  Increasing  wavelength  and  would  be 
even  greater  for  C  ase  1 1  water. 

Variability  In  moorIng-derIved  L„(A)  can  also  be  at¬ 
tributed  to  radiometer  wavelength  shifts  (up  to  ±5  nm 
caused  by  temperature  variations).  Wavelength  shifts  of 
-Fl  to  -f5  nm  by  steps  of  1-nm  were  applied  to  measured 
L  JA,  5m)  and  L^j(A,  15m).  Kl(A,  10m)  was  then  derived 
following  Eqs.  (5a)  and  (5b),  and  L  JA,  Oi  and  L„(A)  were 
computed  using  Eqs.  (6a)  and  (6b),  extraolating  from 
L^j(A,  5m).  Percentage  differences  were  calculated  using 
Eq.  (8),  with  L„(A)2  =  wavelength-shifted  values  and 
=  non-wavelength-shifted  values.  Shifts  of  1-nm 
resulted  In  <2.5%  change  across  all  wavelengths  (aver¬ 
age  difference  was  1.3%),  whereas  -6.5%  changes  oc¬ 
curred  with  a  5-nm  shift  (range  of  2.2-13.2%).  The  blue- 
green  wavelengths  (490  and  510  nm)  were  most  sensi¬ 
tive  to  wavelength  shifts  (data  not  shown).  Results  from 
shifts  of  -1  to  -5  nm  were  the  same  but  of  opposite  sign. 


4.2  Mooring  motion 

Vertical  mooring  heave  effects  on  radiometric  meas¬ 
urements,  which  are  dependent  on  weather  and  sea-state 
conditions,  were  examined  by  deriving  L„(A)  from 
L^j(A,  Oi  extrapolated  using  a  range  of  K(j(A,  z  -F  Az)  and 
Kl(A,  z  -f  Az)  wherez  =  18  m  for  K(j(A,  z)  and  z  =  5,  15, 
and  21  m  for  Kl(A,  z),  and  Az  was  varied  from  0.25  to  3 
m  by  steps  of  0.25  m.  Pressure  sensors  mounted  on  the 
BTM  have  measured  the  maximum  heave  of  the  BTM  as 
about  3  m.  Percentage  differences  were  computed  using 
Eq.  (8),  where  L„(A)2  =  depth-shifted  values  and  L„(A)i 
=  non-depth-shIfted  values.  The  results  show  that  heave 
effects  Increased  with  heave-depth  (Az)  and  depth  of  ver¬ 
tical  attenuation  coefficient,  I.e.  maximum  effects  were 
seen  for  Az  =  3  m  and  for  derivations  using  Kl(A,  18m) 
(Fig.  5).  Fleave  effects  were  also  more  prominent  at  555 
nm  than  other  wavelengths,  although  differences  at  410 
nm  for  Kl(A,  10m)  were  relatively  high.  A  n  upward  moor¬ 
ing  heave  of  1-m  resulted  In  differences  between  depth- 
shifted  and  non-depth  shifted  derivations  of  L„(490  nm) 
ranging  between  2.2%  (3.4%  for  555  nm)  for  Kl(490  nm, 
10  m)  to  3.2%  (7.5%for  555  nm)forKL(490  nm,  18  m). 
Differences  of  6.4%  (9.8%  for  555  nm)  and  9.2%  (20.9% 
for  555  nm)  were  found  for  L„(490  nm)  with  a  heave  of 
3-m  (Fig.  5).  FI eave effects  can  becorrected  by  using  pres¬ 
sure  sensors  mounted  at  the  same  depth  as  radiometers. 

M  ooring  horizontal  motions  related  to  the  watch-cir- 
cle  (I.e.  maximum  distance  a  mooring  travels  around  Its 
anchor)  have  not  been  analyzed  here,  but  horizontal  vari¬ 
ability  can  affect  satelllte-to-mooring  match-ups  (Chang 
and  Gould,  2006).  M  ooring  watch-circles  are  dependent 
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Fig.  6.  Results  from  environmental  variability  analysis.  Simulated  L„(A),  18m)  and  Kl(A,  13m)  computed  using  input 

variables  described  in  Discussion:  Environmental  variability  section  and:  (a)-(c)  Variable  solar  zenith  angles  for  1000  to 
1200  AST,  by  increments  of  15- min  (symbols  denote  every  30- min)  and  (d)-(f)  variable  wind  speeds. 


on  water  depth,  wind  speed  and  direction,  currents,  sea- 
state,  and  type  of  mooring  design  (i.e.  semi-taught  ver¬ 
sus  catenary).  The  BTM  watch-circlefor  the  semi -taught 
design  is  5  km  diameter  whereas  SeaWiFS  and  (M  ODIS- 
aqua)  spatial  resolution  is  ~1  km.  Therefore,  equivalent 
spatial  coverage  of  SeaWiFS  and  mooring  data  may  not 
be  obtained  when  matching  the  latitude  and  longitude  of 
the  mooring’s  anchor  with  a  satellite  pixel.  This  does  not 
pose  a  serious  problem  under  horizontally  homogeneous 
atmospheric  and  oceanic  conditions,  but  may  greatly  af¬ 
fect  match-ups  in  spatially  variable  environments  (i.e. 
coastal  ocean  as  described  by  Chang  and  Gould,  2006). 
Latitude  and  longitude  obtained  from  a  surface-mounted 
global  positioning  system  (GPS)  can  help  alleviate  issues 
associated  with  mooring  horizontal  motions. 

4.3  Environmental  variability 

Mismatches  between  satellite- measured  and  moor¬ 
ing-derived  L„(A)  can  also  be  attributed  to  environmen¬ 
tal  conditions  such  as  solar  angle  variability  due  to  mis¬ 
matched  satel  I  i  te  overpass  ti  mes  and  vari  abl  e  wi  nd  speeds, 
which  are  related  to  transmission  across  the  air-sea  inter¬ 
face  (Eq.  (6b)).  These  factors  were  evaluated  using  the 
radiative  transfer  model.  Hydrolight.  Briefly,  the  well- 
documented  Hydrolight  model  solves  radiative  transfer 
equations  in  water,  based  on  invariant  imbedding  theory 
(M  obley,  1994).  Hydrolight  output  data  relevant  to  this 


study  include  K(j(2,,  z),  Kl(A,  z)  at  z  =  10,  13,  and  18  m, 
and  L„(A).  Necessary  inputs  to  Hydrolight  include  the 
boundary  conditions  and  the  inherent  optical  properties 
(lOPs).  The  lOPs  were  obtained  from  a  historical  Case  I 
water  model,  assuming  two  components:  pure  water  and 
chlorophyll-bearing  particles  with  co-varying  colored 
dissolved  organic  matter  and  detritus.  Pure  water  coeffi¬ 
cients  were  taken  from  Pope  and  Fry  (1997).  Chlorophyll 
concentration  (Chi)  was  assumed  to  be  constant  with 
depth;  Chi  =  0.515  lug  l“\  which  was  the  average  Chi 
measured  by  the  BTM  over  the  period  of  interest. 
Hydrolight  input  boundary  conditions  include  wind  speed, 
solar  angle,  cloud  cover,  downwelling  sky  irradiance,  and 
ocean  bottom  type.  Cloud  cover  was  held  at  0%, 
downwelling  sky  irradiance  was  computed  directly  from 
the  RADTRAN  model  and  waters  were  assumed  to  be 
optically  deep  (i.e.  no  influences  from  the  ocean  bottom). 
Two  different  simulations  were  run  with  these  inputs: 

(1)  Wind  speed  was  held  constant  at  5  m  s^^  while 
solar  angle  was  varied  from  38.0  to  25.7°.  The  angles  were 
obtained  by  calculation  of  solarzenith  anglebetween  1000 
and  1200  AST  by  steps  of  15- min.  In  addition  to  the  sam¬ 
pling  time  of  day,  these  computations  required  latitude 
and  J  ulian  day  (J  D;  equations  (2-2)  and  (2-3)  from  K irk, 
1994).  The  latitude  of  the  BTM  is31.7311°N  andJD  was 
set  at  the  approximate  midpoint  of  the  time  series  of  in¬ 
terest,  JD  =  250  (7  September  1999).  Differences  were 
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evaluated  between  each  solar  angle  and  25.7°  (noontime 
solarzenith  angle),  e.g.,  25.7  compared  to  38.0°,  25.7  and 
35.5°,  etc.  Solar  angle  had  a  great  Impact  on  computa¬ 
tions  of  L„(2,),  as  Indicated  above.  A  2-hr  (1-hr)  differ¬ 
ence  In  time  resulted  In  -15%  (-4%)  differences  In  de¬ 
rived  L„(2,).  Lower  wavelengths  were  slightly  more  af¬ 
fected  by  Increasing  solar  angles  (Fig.  6).  Values  of 
K(j(A,  z)  and  Kl(2,,  z)  remained  relatively  unchanged  (dif¬ 
ferences  of  <1%)  until  solar  angle  Increased  to  more  than 
7.5°  from  solar  noon  or  >1.5  hours  from  1200AST.  Dif¬ 
ferences  then  Increased  to  5-6%.  These  errors  associated 
with  solarzenith  angle  variability  havebeen  reduced  sig¬ 
nificantly  with  the  computation  of 

(2)  Wind  speed  was  varied  from  0  to  15  m  s^^  by 
steps  of  3  m  s^^  while  solar  angle  was  held  constant  at 
30°.  All  other  Inputs  were  the  same  as  for  simulation  (1). 
Differences  reported  were  calculated  between  0  m  s^^  and 
each  Increment  of  Increasing  wind  speed,  e.g.,  0  and  3 
m  s^\  0  and  9  m  s“\  etc.  It  was  found  that  wind  effects 
were  minimal,  with  less  than  3.5%  difference  between 
L„(A)  computed  for  0  m  s^^  and  15  m  s^^  wind  speeds. 
Effects  decreased  with  Increasing  wavelengths  (Fig.  6). 
The  Influence  of  wind  speeds  on  K(j(A,  z)  and  Kl(A,  z) 
were  similarly  small  (<2.2%)  but  Increased  with  Increas¬ 
ing  wavelength  (Fig.  6). 

5.  Recommendations  and  Conclusions 

Our  recommendations  for  the  use  of  moorlngs-of- 
opportunlty  for  ocean  color  satellite  validation 
(ollgotrophic  waters)  are  as  follows: 

•  Comparisons  between  the  pre-  and  post-radlom- 
etercallbratlonsare essential  (M  uellerandAustIn,  2003). 
It  was  confirmed  that  stringent  comparisons  between  the 
two  are  necessary  to  verify  how  radiometer  diodes  might 
drift  during  deployment. 

•  Proper  antl-blofouling  precautions  and  proce¬ 
dures  are  essential  for  successful  data  collection.  In  ad¬ 
dition  to  shutters,  wipers  and  antl-fouling  material,  we 
recommend  that  time- series  data  should  be  Inspected  for 
potential  fouling  signals. 

•  QA/QC  methods,  e.g.,  clear-sky  model  and  pure- 
water  comparisons,  must  be  employed  and  their  usability 
has  been  confirmed. 

•  Although  Instrument  self-shading  effects  are 
relatively  small  In  clear,  oceanic  waters,  correction  fac¬ 
tors  are  relatively  simple  to  compute  and  should  be  used 
as  confirmed  In  the  present  study. 

•  One-hour  averages  of  high  resolution  time  se¬ 
ries  measurements  of  E(j(A,  z)  and  L  JA,  z)  centered  near¬ 
est  to  the  satellite  overpass  time  showed  the  best  results 
compared  to  other  time-averaging  methods,  but  4-hr  av¬ 
erages  around  dally  noon  are  also  suitable.  Averaging 
helps  alleviate  problems  associated  with  wave  focusing. 
Intermittent  cloud  cover  and  other  environmental  vari¬ 


ability  at  site. 

•  Extrapolations  of  LJA,  z)  to  L„(2,)  are  more 
accurate  when  Kl(2,,  z)  Is  used,  as  confirmed  In  the  present 
study.  However,  K(j(A,  z)  can  be  a  reasonable  alternative 
extrapolation  variable.  Here  It  Is  Important  to  clarify  that 
although  K(j(A,  z)  and  Kl(A,  z)  are  completely  different 
optical  properties.  In  some  Instances  when  Instruments 
fall  and/or  only  Instrument  deployments  are  possible  an 
alternative  option  may  be  utilized. 

•  Pressure  sensors  should  be  co-located  with  ra¬ 
diometers  to  enable  corrections  due  to  mooring  heave  and 
Instrumentation  effects,  I.e.  radiometer  wavelength  shifts 
must  be  considered. 

•  A  surface  buoy-mounted  GPS  Is  necessary  for 
accurate  mooring-to-satelllte  pixel  match-ups  In  horizon¬ 
tally  variable  environments. 

•  Radiometers  should  be  placed  within  one  opti¬ 
cal  depth  of  the  sea  surface,  above  the  subsurface  chloro¬ 
phyll  maximum,  but  also  at  a  depth  below  the  Influence 
of  buoy  shadowing  effects.  Procedures  to  address  poten¬ 
tial  buoy  shadow  effects  are  necessary. 

•  At  least  two  upwelling  radiance  and  two 
downwelling  Irradlance  sensors,  I.e.  sensor  redundancy, 
are  recommended  In  the  case  of  sensor  failure. 

While  dedicated  calibration  buoys  such  as  the 
M  OBY,  Y  BOM  and  BOUSSOLE  systems  contribute  to 
direct  satellite  vicarious  calibration  efforts,  we  show  that 
moorings-of-opportunity  can  be  cost-effective  methods 
for  validating  data  collected  by  ocean  color  satellites  when 
using  careful  quality  control  procedures.  Time-averaged 
percentage  differences  between  SeaW I FS- measured  and 
BTM  moorIng-derIved  L„(A)  were  about  20%  at  the  blue 
wavelengths,  decreasing  to  as  low  as  11%  In  the  blue- 
green  to  green  wavebands.  These  results,  however,  do  not 
meet  satellite  ocean  color  goals  of  estimating  L„n(2,)  to 
within  5%.  Nonetheless,  these  relatively  small  differences 
are  quite  remarkable  given  the  extreme  difficulties  In  ac¬ 
curately  measuring  L„(A)  in  situ  and  remotely. 

Our  study  shows  that  existing  and  future  buoy-moor- 
Ing  systems  are  and  will  be  essential  to  ocean  color  valida¬ 
tion  and  vicarious  calibration.  However,  buoy-mooring 
systems  alone  are  not  the  complete  answer  as  they  only 
provide  documentation  of  discrepancies,  not  solutions. 
Future  observational  platform  systems  must  Include,  for 
example,  robust  measurements  relating  to  atmospheric 
optical  and  environmental  conditions.  Optical  data  from 
a  wide  variety  of  oceanic  environments  and  complemen¬ 
tary  Interdisciplinary  data  for  comparative  analyses  are 
essential  to  continuing  Improvements  to  the  quality  of 
satellite  ocean  color  data  and  associated  algorithms  for 
the  derivation  of  environmentally  Important 
blogeochemical  parameters  (lOCCG,  1998,  1999,  2000, 
2004,  2006). 
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